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Abstract 



This article shows the interfacial relation in electrodynamics shall 
be corrected in discrete grid form which can be seen as certain numer- 
ical dispersion beyond the usual bulk type. Furthermore we construct 
a lossy conductor model to illustrate how to simulate more general 
materials other than traditional PEC or simple dielectrics, by a new 
conformal FDTD method which main considers the effects of penetra- 
tive depth and the distribution of free bulk electric charge and current. 
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1 Introduction 



Bulk dynamics and boundary conditions are two essential ingredients 
for solving a concrete physical problem. In computational electro- 
magnetic, it is well known finite difference form of wave equation in 
bulk leads to numerical dispersion which is measured by the factor 
^i^, where x is proportional to the grid size. If the grid length is 
smaller than then the dispersive factor approximates to 1 which 
means the error due to this numerical dispersion can always be ig- 
nored in problem. However, beyond this traditional bulk numerical 
dispersion (BND), we find discrete process in simulation as in Finite- 
Difference Time-Domain (FDTD) approach also generates certain un- 
physical anti-dispersion effects on the boundary which can be defined 
as the interfacial numerical dispersion (IND). If we ignore this effect, 
we shall get the staircase approximation grid representation of the 
curved boundary, which is just one of the shortcomings of traditional 
FDTD method. Instead, we find interfacial relation in usual electrody- 
namics must be corrected in grid variable form due to IND, and some 
treatments in conformal FDTD method before as in perfect electric 
conductor (PEC) or dielectric materials can be seen as special cases of 
the new relation. 

Traditional conformal FDTD (CFDTD) often deal with PEC [1- 
10] and simple dielectrics [11-15]. In former case, the electric field in 
material can be ignored due to the penetrative depth is far less than 
one grid length due to the lossy parameter p ~ ^ 1 and in later 
case simple denote near lossless, i.e. p = ^ <C 1. 

For general dispersive and lossy media like Lorentz model, the re- 
lated Conformal FDTD approach is difficult due to the permittivities 
in the interfacial relation involve complex dependence of the frequency 
and in corresponding time-domain algorithm this shall turns to be in- 
volve with high order time derivatives. 

However, for a large class of lossy dielectrics, whose penetrative 
depth i5 is between PEC and perfect electric insulator (PEI), more ex- 
actly in which the lossy parameter p ~ ^ \ oi the 5 is about several 
grids' length, if we can ignore the influence of the resonance frequency 
in the problem or the involved electromagnetical frequency is far away 
from the material resonance frequencies, we can use certain conductiv- 
ity model as the dispersive model and get a (theoretical) more effective 
CFDTD method called lossy conductor CFDTD (LCCFDTD), which 
consider the effects of IND, surface/volume free electric current /charge, 
penetrative depth and conformal angle (defined as the angle between 
the cross-section of interface and the grid line) in more physical way. 
Especially our method is more suitable for the FDTD simulation of 
the problems like the design of microwave dark chamber with strong 
absorber materials. 
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2 Grid representation of interfacial rela- 
tion 



Consider the 2-dimensional TE case below without loss of generality. 
First we define the angle 9 between the tangent line at the interfa- 
cial point and the grid line as the conformal angle. The grid line at 
the boundary is divided into two parts on which two electric fields is 
uniform distributed as E in air and E' in materia. General case of 
arbitrary two medias is similar. 

For PEC model, the inside field can be assumed to be vanished. But 
for general materials or even as in nondispersive, lossless dielectrics, to 
obtain a conformal algorithm instead of naive staircase approximation, 
we have to discuss the interfacial relation between E and E' . In con- 
tinuous electrodynamics theory we have following boundary relations: 



where a is line density of free surface electrical current and p is surface 
density of free electrical charge. Here we only consider the simplest case 
in which a, p vanish and permittivity is a real const. In next section 
we may show how to deal with nonzero a, p and complex permittivity 
for a more general model of materials. 

As we know, for certain frequency only two of above interfacial 
relations are independent and here we use (1) (3), i.e. 



This interfacial relation indicate E, E' can't be along same direction 
through the interface, i.e. the refraction phenomena exist or for general 
case if the permittivity as function of the frequency, electromagnetical 
waves of different frequency shall refract to different directions, i.e. 
dispersive phenomena appears. However, in discrete grid simulation, 
E, E' have to be along same grid line and this in fact introduces 
extra error. Like the discrete treatment for wave equation in bulk 
will produces numerical dispersion, the effect of this elimination of 
physical refraction or dispersion on interface due ro simulation method 
as FDTD can be seen as certain interfacial numerical dispersion (IND). 
Hence we should take one of two fields on the same boundary grid line 
as projection of the refraction field along the grid line, take E for 
instance and the relation between E, E' give the grid representation 
of physical interfacial relation (1), (3). 



K 



(1) 

(2) 
(3) 
(4) 



E,, = E'^, E±_ = ErE'j_, Sj- = e/eq. 
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By a simple geometrical analysis indicated by figure 1 above, we 
find this grid interfacial relation, 

E={cos^9 + erSm^e)E' ^eriO)E'. (5) 

It's direct to see several special correct results. First take Sr — I i-e. for 
grid in same materia, E — E' as our basic assumption in simulation. 
Second, consider = 0, i.e. the coincident case of grid choice, E = E' , 
as same as traditional FDTD. Third, if 9 — ^, as the translation 
case for coincident grid choice, we get E — ErE' or equivalently as 
E„=E^,=0, Dj_=D'j^. 



3 Lossy conductor model and conformal 
FDTD 

To deal with more general materials as needed in most applications, 
we construct a lossy conductor model and for simplicity, we assume 
the boundary of materia is flat, i.e the tangent line at the interfacial 
point is coincident with the boundary. New model main considers 
a reasonable generalization of the surface distribution of free electric 
charge and current in well conductor model (WEC) via a special dis- 
crete treatment. 

First we introduce a kind of lossy conductor conformal grid (LCCG), 
and take a 2-dimensional cross-section of a rotated planar interface for 
example as figure 2 shows. Except usual field nodes, this grid have 
four material points a, 6, c, d and a new grid parameter B called grid 
penetrative number. For the grid outside of the material B — and 
on the geometrical boundary of media B = 1. For the grids in the 
material, B — mm{[Di/ds]) + 1, where is the distance between the 
grid's center node and the tangent line of interface, [x] is the maximal 
integer less than x, ds denote the size of grid. In fact we can take its 
relative coordinate of the nearest interfacial grid as its B value. 

Only two material points in a grid are independent which are de- 
termined by the dividing character of the boundary grids, while other 
two points if needed can be obtained via translation of the line which 
connects the two independent points as figure 2 shows. In fact like 
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traditional subgrid method which introduces four field nodes in a grid 
to construct the relative small grids, our new conformal method also 
consider new sub nodes which describe the geometry of material bound- 
ary and this process don't introduces new field nodes to increase the 
complexity of computation. 

The start point of conformal method is the integrated form of 
Maxwell equation below: 



dt 



[ iiH-dS = - <p E-dl. (6) 

Js JdS 



For the case of figure 2, by discretion of above equation we have iterated 
equation for magnetic field H, 

H^HiJ) =HrHi,j) + J^^^iKiiJ + + E:{i,j + \)Ui,j) 

- E:{t,j - \)lAhJ - 1) - E'-{t,j 1) 
+ E^ii - l,j)ly{i - + E^-{i - - 1,3) 

- E^{i + \j)ly{i,j) - E'^^H + \,j)l'y{i,3)] 

Note we can just take some terms in above formula to be zero for 
other cases of division. In the situation of figure 2, we also have some 
constraints on the length variables due to two independent materia 
points in a grid, 

l'y{i,j) = ly{i - UiJ - 1) = l.ihj) ■ i^l^ZYJy 

I +1' =Ax I +1' =Av Mm) 

+ Ax,ly + l^ l'AiJ)~l'yii,jy 

It's easy to sec that the usual PEC conformal algorithm is just as the 
special case for which S{i,j) denote the area of the outside cell and 

E.{i,j + l)= E',{i,j + 1), Ey{i - 1, j) = E'yii - 1, j); 

E',{i,j-^) = 0, E'y{z+^,j) = 0. 

By the symmetry of electric grids and magnetic grids in FDTD 
method, it's not necessary to do conformal treatment for electric grids, 
i.e. the iteration of electrical field is same as in traditional FDTD 
which is just by discretion of another Maxwell integrated equation: 



as 



H-dl = ^^ I D-ds+ I J-ds. 



Note in our lossy conductor model, the conductivity a is not zero which 
require the term of free electric current. 

To deal with the effects of the free electric charge and penetrative 
depth, we take a special discrete treatment for the distribution of free 
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bulk electric current in our grid system, which assumes the free electric 
charge only distribute uniformly on the transverse grid line It and the 
electric current just conduct along the perpendicular direction of the 
grid plane. Thus no electric charge and current on the grid line in 
our model and especially the transverse line It in one grid just denote 
the unit transverse line which appears in the definition of free surface 
electric current a in (2). 

By a simple deduction like for the lossless model, we can get the 
generalization of formula (5) in lossy conductor model, 

E = (cos2 9+{er- — + sin^ e)E' 

= (cos^ + Ere sin^ 6)E' = e„(6»)£'. 

where |a| = gE'It as the electric current through the unit transverse 
line It and v is the velocity of the electromagnetical wave in the ma- 
terial and clearly v = ^^^^ = where the refraction lossy 

parameter L = \/^^^^^, p = -f^, Or = cr/eo- Finally we obtain 
the grid interfacial relation of lossy conductor model: 

E = (cos^ e + {Sr-^ + ^1^1^^^) sin2 e)E' (7) 

To apply above formula in conformal FDTD scheme, we have to 
make some approximations and furthermore turn it to the iteration 
relation in time domain. 

• Half order and IND corrected approximation (HIND): It — 
(ignore the distribution of free electrical current) 

E={cos^e+{er-^)sm'^0)E' (8) 

dE dE' 

— = (cos^ e + Er sin^ 6*) — (Tr sin^ e-E'. (9) 

at at 

• First order and IND corrected approximation (FIND): 

It^'O, Lc^l 

E = (cos^ e + {Er- — + ^1^1^^) sin2 e)E' (10) 
juj c sm 9 

dE , o , iTCfrx/ErUr s 9 ^dE' 9 , , , 

ot csmO at 

• Second order and IND corrected approximation (SIND): 
It^G, L~1 + ^ 

E = (cos^ e + {Er-^ + ^^^''^ (1 + ^)) sin^ e)E' (12) 

= cos2 e + {£r + ^-^) sm' ^ W - sm' 9 ■ — 

at-' csm9 at'' at 

-'^..oe: (.3) 
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• Geometric boundary conformal FDTD: 

Only consider above conformal treatment for the geometrical bound- 
ary grid in which B = 1. Thus in the iteration equation of magnetic 
field we have 

EA^.J - ^) = Ki^,J - I), Eyii+^,j) = E'y{i+^,j) 

a.nd Ex{i,j + ^), Ey{i-~ ca,n he dctcimmcdhy E'^{i, j + E'y{i — 
i, j) via the discrete form of the equation (9), (11) and (13). Fori? > 1, 
we set E = E' which means we have ignored the effect of penetrative 
depth in this scheme. 

• Penetrative boundary conformal FDTD: 

Consider above conformal treatment for the penetrative boundary 
grid in which 1 < B < [^], M = min(Aa;, Ay). For such grids we 
get all E from E' via discrete form of (9), (11) and (13). The grid 
with i? = [^] + 1 will be treated as geometric boundary grid while for 
B > [^] + 1 we just set E = E'. Especially for PEC case, [^] = 
and = as in traditional way. 

Although the penetrative depth S is an important physical quantity 
which can be determined by the electromagnetic parameters of the 
material exactly, we can take it as an adjustable parameter called lossy 
depth Si for practical use. It can be set by the precision and cost of 
computation in a problem and penetrative depth i5 could just be its 
default value. 

4 Discussion 

We know electrodynamics manifest Lorentz symmetry or Maxwell equa- 
tions in continuous spacetime is covariant under the translation and 
rotation, while in discrete grid spacetime such Poincare symmetry is 
broken. For flat boundary material, we find formula (5) considers the 
exact effect of conformal angle in the interfacial relation such that the 
simulation physical result is independent of grid system in the sense 
of global translation and rotation, or in other words for flat boundary 
case the shape of material in simulation is invariant of grid translation 
and rotation which is more physical unlike in traditional FDTD. In this 
manner we could say our IND corrected conformal FDTD have restored 
certain 3-dimensional space symmetry and especially it manifests the 
consistence between the conformal and nonconformal algorithm. 

There are several directions to generalize above section's result. 

First, if the boundary is curved, the tangent line at the interfacial 
point shall not be along It- One simple way to deal with this case is 
to define the angle between the tangent line and transversal line It as 
relative conformal angle 6r and the conformal angle still to be the one 
between It and the grid line, li 9r < |, such boundary's treatment 
can be approximated by transversal line boundary as done in previous, 
while 9r > we'd better approximate the boundary in this direction 
by the complete grid line, i.e. E = E' . 
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Second, our LCCFDTD method only considers the electrical loss, 
it seems the magnetic loss and some special anisotropy or nonlinear or 
dispersive medias also can be treated in similar way [16]. It may also 
have applications in other numerical computational areas. 

Finally, we may consider the stability problem as in traditional con- 
formal FDTD method [17-18] and an improved absorb boundary may 
be obtained by the combination of a large lossy conductor material and 
the ideas of PML. Of course all above discussions are main in theoreti- 
cal manner and it can be verified by some numerical experiments. The 
results may be similar with [14], i.e. the simulation of rotated problem 
should has good agreement with the coincidence problem as reference 
solution. 
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